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Abstract 



In the light of the recent WMAP results we update the constraints on a 
class of non standard BBN models with a simultaneous combination of non 
standard neutrino distributions and extra effective number of neutrinos in 
the expansion rate. These models can be described in terms of the two pa- 
rameters AA'^J"*, constrained by the primordial Helium abundance Yp mea- 
surement, and ANu, constrained by a combination of CMB and primordial 
Deuterium data. Small deviations from standard Big Bang Nucleosynthesis 
are suggested. Different non standard scenarios can be distinguished by a 
measurement of the difference AnI" = ANl°^ — AN^. From the current 
data we estimate AnI" ~ —1.4^5*45 mildly disfavouring solutions with a 
low expansion rate, characterized by AnI" = and negative AN^. Active- 
sterile neutrino mixing could be a viable explanation only for high values 
of Ij, > 0.24. The existence of large positive neutrino chemical potentials 
~ 0.05, implying AA''^ ~ 0, would be a possible explanation of the data 
within the analyzed class of non standard BBN models. Interestingly it would 
also provide a way to evade the cosmological bounds for 'class A 3+1' four 
neutrino mixing models. A scenario with a decaying sterile neutrino is also 
considered. 
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I. INTRODUCTION 



In a previous paper [1] (see also [2]) we showed how the new CMB measurements of 
the baryon to photon ratio, rj, are able to put stringent constraints on a large class of non 
standard BBN models where, together with a usual variation of the expansion rate due to 
the presence of extra degrees of freedom, distortions of the electron neutrino distribution 
are also present. This class of models can be described in terms of two parameters [3]. A 
first one is the usual extra cfi'cctive number of neutrinos, modifying the standard expansion 
rate, AA^^ = (Px / Po) ~ 3] , where px is the energy density of the X-particle species, 
including the three ordinary neutrinos plus possible new ones, and po — (77r^/120) is the 
energy density of one standard neutrino species. The second one is the total extra effective 
number of neutrinos AA^^°* defined, in terms of the primordial ^He abundance Yp, as: 

AAr*°t = [Yf^'^iv, AN^, 5UJ - Yp'^^""] /0.0137 . (1) 

The difference ANjf^ — ANj^ is a quantity that, in the class of models that we are consider- 
ing, has to be entirely ascribed to the effect of deviations of the electron electron neutrino 
distribution from the standard Fermi-Dirac with zero chemical potential, 5f^^ — f^^ — f^^. 
If Sf^^ = then AA^*°* = AA^^ and simply [4]: 

r/^^(r;, AAT^^ 6f.^ = 0) ^ Y^'^^'^iv) + ^{rj) AN^ (2) 

with 7] the baryon to photon ratio in units of 10~^°. Using the expansion given in [4], we 
calculated that 7(77) ~ 0.0137 over the pertinent range 7] = 3.5 — 10. The standard BBN 
prediction for Yp is described by the following expansion around rj = 5 [4]: 

ySBBN _ Q 2456 + 0.01 In (77/5) (3) 

The presence of non zero 6f^^ affects mainly Yp, while its effect can be safely neglected 
in the Deuterium abundance, {D/H), also considering that we will be interested in small 
deviations. With this approximation the D/H abundance is described by the expression [1]: 

(l + aAA^^)2 (4) 

with /3 ~ 1.6 and a — rl^/ Ag^Q^'^) ~ 0.135, where r,^o and are, respectively, the 
standard neutrino to photon temperature ratio and the number of degrees of freedom at 
present. With these expressions a simultaneous measurement of (D/H), Yp and 77 can be 
easily translated into a 'measurement' of AA^*°* and AA^^. We used in [1] both high [5] ^ 

y^f^P = 0.244 ± 0.002 (5) 

and low values Fp'^^P = 0.234 ±0.003, while we used (D/iJ)^^P = (3.0 ±0.4) x 10"^ [7]. For rj 
we used the DASI and BOOMerANG result [8] rj^^^ = G.O^oi- From low values of Helium 



{D/Hy'"\ri,ANI^) ~ 3.6 • 10"^ (77/5) 



We indicate 68% c.l. errors for all quantities unless differently indicated. 
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and assuming gaussian errors, we obtained AiV*°* = —1.05 ± 0.25, while from high values 
of Helium we obtained ANl°^ = —0.3 ± 0.2. Using the primordial Deuterium abundance 
measurement, from the expression (4), we could estimate AA'^^ obtaining AiV^ = 1±4. These 
results were implying, at 3 cr the bounds [1] AA^*°* < 0.3 and AA^^ < 13. In particular the 
bound on AA^*°* was used to conclude that, for neghgible neutrino asymmetries, all four 
neutrino mixing models are in disagreement with cosmology and thus ruled out. This has 
been then also confirmed by the improved solar and atmospheric neutrino data from the 
SNO [9] and SuperK [10] experiments [11]. In the next section we will update these results 
in light, mainly, of the recent results from the WMAP experiment [12] and we will see how 
the data suggest possible deviations from a standard picture. 



II. UPDATED REFERENCE VALUES AND RESULTS 

The WMAP collaboration finds fib = 0.0224 ± 0.0009 [12] corresponding to: 

r^c^AfB ^ g 15 ^ 0.25 (6) 

This measurement is so precise that now, when estimating ANl°^, the experimental error 
on Yp is dominant compared to that one on rj. Using high values of Y^^^ we find at 1 a: 
AA^*°* _ _o.35 ± 0.15. This means that now a 3a range is given by 

-0.8 < A7V^°' < 0.1, (7) 

implying a quite more stringent upper bound compared to the pre- WMAP value. Even 
using the range of values 

y;^P = 0.238 ± 0.002 ± 0.005 , (8) 

that is a compromise between low and high values and takes into account the discrepancy 
as a sistematic uncertainty [13], we find 

ANl"' = -0.8 ± 0.4 , (9) 

implying a 3 cr range 

-2.0 < AATtot ^ Q 4 ^10) 

Both results confirm our previous conclusion for which AA^*°* as high as 1 is highly dis- 
favoured, thus ruling out all four neutrino mixing models in the case of negligible neutrino 
asymmetries [1]. However now both results seem to point out, at 2(7, to a negative value of 
AiV*°*, suggesting the presence of non standard BBN effects. We can also update the esti- 
mation of AA^^ using the new r) measurement from CMB and a new primordial Deuterium 
abundance measurement [14], {D/HY^'p = (2.781°;^^) x 10-^ finding {AN^)^^^ = 0.7±2.1. 
As already anticipated in [1], the error has been highly reduced by the great improvement 
in the rj determination from CMB and it is now dominated by the error on D/H. However, 
differently from the determination of ANjf^ from Y^^^, better future determinations of rj 
(for example from new WMAP data or from Planck) can still further reduce the error on 
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(AN^)^^^ from the current 2.1 down to 1.5. It is interesting to note that the value from 
BBN is comparable to the direct determination from CMB. In [15], combining the WMAP 
data with the 2dF redshift survey and using the value on the Hubble constant from the 
HST Key Project, h = 0.72 ± 0.08 [16], the authors find {ANPf^^ = O.SlJj. Assuming 
that, between the nucleosynthesis and the recombination time, the quantity AA^^ does not 
change ^ and thus that (AN^)^^^ — [AN^)^^^ ^ one can then combine the two values. 
We will still assume gaussian errors for a qualitative estimation ^ and in this way we find a 
CMB-Deuterium combined value 

In this way we get a much more stringent 2 a (3(t) range: 

-(1.8) 1.0 < (AAr^)^^^+^/^ < 3.4 (4.8) . (12) 



III. POSSIBLE SCENARIOS 

These new results show that deviations from Standard BBN, if they exist, are small. This 
means that Standard BBN is in any case, in first approximation, a very good description 
of all data. This result is mainly due to the fact that the Deuterium abundance is in very 
good agreement with the CMB prediction. At the same time the measured primordial 
Helium abundance, 1^, suggests the possible presence of small deviations whose detection 
is now possible mainly to the great precision of CMB in measuring the baryon asymmetry. 
However, for an assessment of such a hint, it will be necessary to reduce the large sistematic 
uncertainties on Yp and it will be also necessary to investigate even more accurately on the 
robustness of the rj determination from CMB. In the following we will assume that such a 
hint is suggestive of non standard BBN effects and we will discuss some possible scenarios 
that could explain these deviations. An important role in our discussion is given by the 
quantity ANl^ = AA^*°* - AA^^. Prom (11) and (9) we can estimate: 

AN I'' ~ -l-4l?i (13) 



A. Low expansion rate 

A minimal possible way to interpret the data is to assume that there is no efi^ect due to 
electron neutrino distribution distortions and thus AN^" = or equivalently AA^^°* = AN^. 
In this case one can combine the result (9) from and the result (11) from Deuterium plus 



^See [1] and [17] for discussions and examples in which (AN^)'^^^ ^ (AN^)^^^ . 

^From the Ukelihood distribution given in [15], this does not seem to be a very good approximation 
at values larger than the central one, while it is reasonably good for smaller values. 
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CMB, getting (AiVP)CMB+D/ff+yp ^ _q_q+oao_ rpj^jg ^.^g^^^ ^^^j^ suggest a negative value 
of AN{^, mainly due to the low value of Yp, implying a highly non standard modification of 
the expansion rate during the BBN time, more precisely a lower expansion rate. Usually the 
presence of new particle species would lead to a higher expansion rate and therefore such 
a possibility should rely on some drastic change of the radiation dominated picture during 
the BBN period. However note that, from the Eq. (13), the measurements mildly favour 
a value ^Nl" ^ and so this scenario is mildly disfavoured from the data (at almost 90% 
c.L). 

B. Degenerate BBN 

A well known modification of the Standard BBN is to introduce neutrino chemical poten- 
tials in the thermal distributions [20], corresponding to have pre-existing neutrino asymme- 
tries or generated at temperatures T > 10 McV by some unspecified mechanism. An electron 
neutrino chemical potential (^e = l^e/T) would yield AA^*°* ~ — 16^e- The observed Yp (cf. 
(8)) would then be explained by having 

= 0.05 ± 0.025 . (14) 

It has been shown in [18], extending the results of [19], that the existing information on 
neutrino mixing makes possible to conclude that before the onset of BBN arbitrary initial 
neutrino chemical potentials would be almost equilibrated in a way that ~ ~ ^g. The 
presence of chemical potentials would thus correspond to 

~ 3 X 10"^ < AAr^°* (15) 

Therefore in this scenario the expansion rate would be practically standard and the devia- 
tions would entirely arise from non standard electron neutrino distribution. 

C. Active sterile neutrino oscillations 

Let us assume now that at temperature T ^ lOMeV all neutrino asymmetries are 
negligible, for example of the order of the baryon asymmetry. It has been shown in many 
papers [21] that a small mixing betwen active neutrinos and new light sterile neutrinos can 
generate ordinary neutrino asymmetries and thus negative values of AA/"/" together with 
AA^^ > 0. In a simplified two neutrino mixing the value of AN^" is highly dependent on the 
value of the parameter Am^^ = — m^. Usually the possibility to introduce active-sterile 
neutrino oscillations was motivated by the LSND anomaly [23] . However an explanation of 
the LSND anomaly in terms of active-sterile neutrino oscillations, compatible with the solar 
and atmospheric neutrino data would yield, as already mentioned, ANjf^ = ANj^ ~ 1 [1] (see 
also [24]). At the same time the new WMAP bound on the neutrino masses, < 0.23 eV 
[12], is now also incompatible with such an explanation of the LSND anomaly [25], except 
for one constrained exception [26] . The possibihty to generate a negative AA^*°* requires a 
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negative value of Am^ = ml — mf and very small mixing angles (sin^ 2^ -C 10 [21,1]). 
Values of < 0.23 eV imply thus lAm^J < 5 x lO^^eV^ In [22] it was shown how such 
maximum value, together with very small mixing angles, would produce AN^" > —0.3. 
For an inverted full hierarchical case the corresponding \Am'^g\ ~ 10~^eV^ and in this case 
AA^/" ~ —0.13. These values have to be considered as maximal because in the reality one 
should consider a full multi-flavour mixing and, though full calculations are still missing, one 
can expect that part of electron neutrino asymmetry is actually shared with the other two 
flavours. This means that the small effect could reconcile the observed tjb from CMB only 
with high values of Yp (cf. (5)). In a two neutrino mixing small positive values of AA^^ are 
also possible, for larger mixing angles, but this would go at expenses of jAA^/^l, making it 
even smaller [1] . Having more than one sterile neutrino flavour would make possible to have 
ANI" ~ —0.3 and positive AA*"^ but in this case the total AA^*°* would be larger than —0.3. 
This possibility is however interesting, since it would be a way to distinguish active-sterile 
neutrino oscillations from a degenerate BBN scenario. Another way would be the detection 
of the effects of a possible formation of neutrino domains [27] , hke inhomeogeneities in the 
primordial Deuterium abundance [27] that would give rise to gravitational waves [28] . 



D. Degenerate BBN and class A 3-|-l models 



This is an intriguing variation of the pure degenerate BBN scenario. Suppose there 
are both large chemical potentials and also a mixing of new sterile neutrino flavors with the 
ordinary ones. If the chemical potentials are of the order given by the Eq. (14), then, even for 
maximal mixing, the sterile neutrino production prior the onset of BBN would be suppressed 
[29] and consequently the final value of AN^^ would be the same as in the degenerate BBN 
scenario, while AA"^ can be in principle slightly higher because of a initial sterile neutrino 
highly diluted abundance. In this way it is very interesting that, as already noted in [1], 
the cosmological bound on four neutrino mixing models can be evaded. Moreover this same 
conclusion applies also to the (current WMAP or any future one) bound on the sum of 
neutrino masses applied to so called 'class A 3+1' four neutrino mixing models. They are 
such that the highest mass eigenstate is almost coinciding with a new sterile neutrino flavour 
and separated from the three lighter ones, almost coinciding with the ordinary ones (see [1] 
for references and details), by the LSND gap. In this case the sterile neutrino contribution 
to the fraction VLy/VLrn would be negligible and the bound on the sum of neutrino masses 
would apply only to the three active mass eigenstates whose total mass is the same as in 
ordinary three neutrino models. Note also that among four neutrino mixing models these 
are the only ones to be still marginally consistent with neutrino mixing experiments [30]. 



E. Decaying sterile neutrino 

CP symmetrical decays of a sterile neutrino with a mass nist ^ into electron neu- 
trino plus, for example, some unknown scalar, with a lifetime r, could yield a positive AA"^ 
and at the same time a negative AA"*°* analogously to the decaying MeV — Vr mechanism 
of Hansen and Villante [31], but with some important differences. A sterile neutrino abun- 
dance produced before the quark-hadron phase transition, at temperatures T » 100 MeV, is 
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necessarily highly diluted compared to that of ordinary neutrinos. This can be rc-gcncrated 
at a lower temperature ~ 15 McV (mst/cV)^/"', implying rrist ^ 1 KcV, by a possible neutrino 
mixing with ordinary states (see [1] for details and references). If the decaying temperature 
T{t — t) is approximately comprised in a window (0.5 — 5) MeV, elastic scatterings can 
partially or totally kinetically equilibrate the excess of produced electron neutrino and anti- 
neutrinos but its total number cannot be totally destroyed by the partially or completely 
frozen annihilations before the freezing of neutron to proton ratio. This symmetric excess 
of electron neutrinos and antineutrinos would yield a negative AA*"*"*^ that can agree with 
the value (9). Note that since the decaying neutrinos are sterile, their decays would anyway 
occur out-of- equilibrium even though in the ultrarelativistic regime. An interesting possibil- 
ity is the case that the sterile neutrino is the LSND neutrino of class A 3-1-1 models, with 
mst ~ a/ Am|^^^, like for the degenerate BBN scenario in 111 D. Now however AN^ < 1, 
because the decays can only partly destroy what generated by the mixing. This scenario 
clearly suffers from fine tuning between the lifetime of sterile neutrinos and the time window 
between freezing of Vg annihilations and of the neutron to proton ratio. A way to circum- 
vent this problem is to allow decays to be CP asymmetrical in a way to realize a sort of 
Fukugita-Yanagida leptogenesis at low temperatures. In this case it is enough that the life 
time is shorter than the n/p freezing time (~ 10 sec). This case would be very similar to the 
scenario HID but with AN^^ < 1, the exact value depending on the mixing and on the life 
time. Note that hke in HID the cosmological bound on the sum of neutrino masses would 
be also evaded, as pointed out in [32]. 

IV. CONCLUSIONS 

In future years a better understanding of sistematic uncertainties in the measured Yp 
could strengthen or disprove the hint of non standard BBN effects. At the same time im- 
proved data from CMB experiments should both be able to measure AA^^ with a precision 
of ~ 0.1 [33,17] and make even more robust and precise the determination of rjB- If the 
primordial Helium anomaly will be confirmed, implying negative ANjf^ < 0, then a key 
quantity in discriminating among different explanations is the difference AA^*°* — AN^. If 
this will prove to be not zero and negative, then low expansion rate scenarios will be ruled 
out, as already mildly suggested from current data, and a scenario with large chemical po- 
tentials would be a possible explanation if at the same time AiV^ ~ O{10~^ — 10~^) (maybe 
detectable in a very optimistic case [33]). In the case that lAA/"*"*! < 0.3, then active-sterile 
neutrino mixing can be a viable explanation too and if this is also accompanied by a positive 
value of AiV^, then it will be actually favoured, since degenerate BBN would be ruled out. 
We also pointed out that the degenerate BBN scenario could receive support from neutrino 
mixing experiments. This because large neutrino asymmetries would make 'class A 3-1-1' 
four neutrino mixing models a viable cosmological solution without any limitation from the 
bound on the sum of the neutrino masses. In this case one should receive a confirmation 
of LSND from MiniBoone [34] that would realize a nice consistency between current cos- 
mological data (but we need a better understanding of primordial Helium measurements) 
and neutrino mixing experiments. We would be left with the tough theoretical problem 
to understand the origin of large neutrino asymmetries. Current knowledge excludes the 
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nice possibility of active-sterile neutrino oscillations themselves, but maybe a further in- 
vestigation could change such a conclusion, in particular considering that full multi flavor 
active-sterile neutrino mixing calculations are still missing and that the role of phases in 
three neutrino mixing has never been studied [1]. One possibility is that the sterile neutrino 
decays generate the needed electron neutrino distortions. This case could explain the cur- 
rent central values of AiV*°* (cf. (9) ) and AA^^^ (cf. (11)) but only future more accurate 
determination will allow to distinguish among the different scenarios first of all between the 
standard scenario and possible non standard ones. 
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